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Azimuthally-sensitive femtoscopy for heavy-ion collisions at RHIC and LHC is explored within 
the approach consisting of the hydrodynamics of perfect fluid followed by statistical hadroniza- 
tion. It is found that for the RHIC initial conditions the very same framework that reproduces the 
standard soft observables (including the transverse-momentum spectra, the elliptic flow, and the 
azimuthally-averaged HBT radii) leads to a proper description of the azimuthally-sensitive femto- 
scopic observables - we find that the azimuthal variation of the side and out HBT radii is very well 
reproduced for all centralities, while the out-side correlation is somewhat too large for non-central 
events. Concerning the dependence of the femtoscopic parameters on kr we find that it is very 
well reproduced for the out and side radii, and fairly well for the long radius. The model is then 
extrapolated for the LHC energy. We predict the overall moderate growth of the HBT radii and 
the decrease of their azimuthal oscillations. Such effects are naturally caused by longer evolution 
times. In addition, we discuss in detail the space-time patterns of particle emission. We show that 
they are quite complex and argue that the overall shape seen by the femtoscopic methods cannot 
be easily disentangled on the basis of simple-minded arguments. 
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I. INTRODUCTION 



Femtoscopy provides detailed information on the dy- 
namics of systems formed in relativistic heavy-ion colli- 
sions (for a review and literature see [l|T). Inparticular, 
the Hanbury-Brown-Twiss (HBT) [111 H Hf intensity 
interferometry exhibits sensitivity to shape and flow of 
the medium at freeze-out, which in turn reflects the con- 
ditions throughout the evolution, from the formation till 
the cease of interactions. Together with other observed 
quantities, such as the ratios of particle abundances, mo- 
mentum spectra, the elliptic flow coefficient, or other 
correlation data, femtoscopy contributes to the growing 
precise knowledge of the dynamics of the system. Az- 
imuthally sensitive HBT interferometry (azHBT) [S 0] , 
which is the subject of this paper, brings in information 
on the dependence of shape and flow on the azimuthal 
angle <f>. This information is complementary to the data 
on the transverse-momentum elliptic coefficient v-i Q, 
which is a measure of the azimuthal asymmetry of the 
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flow. Thus, a simultaneous description of azHBT and v 2 
verifies consistency of the theoretical description. 

In this work we use what we call the standard approach, 
consisting of ideal-fluid hydrodynamics followed by sta- 
tistical hadronization. Numerous calculations have been 
done in this framework, with the common difficulty [| 
of simultaneously describing femtoscopy and other signa- 
tures of the data. More precisely, the RHIC HBT puz- 
zle [3, IS ES EH refers to problems in reconciling the 
large value of the elliptic flow coefficient, t>2, with the 
HBT interferometry in calculations based on hydrody- 
namics 

US ES El EI Eft 

In the standard approach 
the description of V2 requires longer hydrodynamic evo- 
lution times, while the HBT radii, sensitive to the size 
and lifetime of the system, are properly reproduced when 
the lifetime of the system is short. The two require- 
ments are in conflict. Due to the obstacle of the RHIC 
HBT puzzle, only few calculations or physically moti- 
vated par ametrization of the azHBT quantities have been 
made [ll El ES IM E3, US H • Moreover, in our view 
the comparison to data of the azimuthal dependence of 
the correlation radii makes sense only when the values 
averaged over <f> are also properly reproduced. 

Recently, we have accomplished a successful uniform 
description of soft observables at RHIC, including the 
HBT radii, within the standard approach [24J. The es- 
sential ingredients of this analysis are the Gaussian initial 
condition for hydrodynamics, early start of the evolu- 
tion, the state-of-the art equation of state with smooth 
crossover, and the use of THERMINATOR [H with all res- 



2 



onances from SHARE [26[ incorporated to carry out the 
statistical hadronization at the freeze-out surface of tem- 
perature 145 MeV. The interplay of these elements re- 
sulted in a simultaneous description of the transverse- 
momentum spectra of pions, kaons and protons, the V2, 
and the HBT correlation radii of pions, including the 
basic azHBT signatures. The agreement of this quality 
cannot be achieved with the initial condition obtained 
from the typically used Glauber models. In that case the 
HBT radii, in particular the ratio of out to side radii, is 
reproduced at the level of 20% only [27]]. The reader is 
referred to Ref. [24[ for the details. 

In this paper we present a systematic study of the az- 
imuthally sensitive HBT radii in the model of Ref. [3] , 
including the dependence on centrality and the transverse 
momentum of the pair, fc-r. We find that the azimuthal 
variation of the side and out radii is properly reproduced 
for all centralities, while in the case of the out-side cor- 
relation it is somewhat too large for non-central events. 
Concerning the dependence of the femtoscopic parame- 
ters on kx, it is very well reproduced for the out and side 
radii, and fairly well for the long radius. We stress that 
the present study involves no parametric freedom, as all 
parameters have been fixed in the global fits of Ref. [24| • 

Next, we extrapolate the model to the LHC energy, 
where predictions for the HBT and azHBT quantities is 
made. In addition, we present a detailed analysis of the 
space-time patterns of the pion emission. We argue that 
this emission, consisting of surface and volume parts, ex- 
hibits a rather complex behavior. We argue that the 
determination of the overall shape of the source seen by 
the femtoscopic methods is not straightforward and re- 
quires detailed simulation, such as the one performed in 
this work. 

We use c = 1 throughout the paper. The label RHIC 
denotes the AuAu collisions at the highest RHIC en- 
ergy of ^snn — 200 GeV, while LHC corresponds to 
the PhPb collisions at ^/snn = 5500 GeV. 



II. THE FRAMEWORK 

In this section we describe the essential elements of 
our method to the extent they are necessary for the 
comprehensive presentation of the new results. More 
details concerning the hydrodynamics can be found in 
Refs. 0, [H, [29|, while the method used for femto- 
scopic calculations has been presented in great detail in 
Ref. 11. 



A. Initial condition 

As reported in Ref. [24[ , the use of the initial condition 
for hydrodynamics of the Gaussian form, 
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FIG. 1: (Color online) In-plane and out-of-plane sections of 
the two-dimensional energy-density profiles for c = 30 — 40% 
obtained from the mixed Glauber model described in the text 
(dashed lines) and the Gaussian parametrization fT} used in 
this work (solid lines). Both profiles are normalized to unity, 
and in addition {x 2 ) = a 2 , (y 2 ) — b 2 . 



where x and y denote the transverse coordinates, leads to 
a much better uniform description of the data for the px- 
spectra, V2, and the pionic HBT radii compared to the 
use of the standard initial condition from the Glauber 
model. 

The width parameters a and b depend on centrality. In 
order to use realistic values we run the GLISSANDO [3l[ 
Glauber Monte Carlo simulations which include the ec- 
centricity fluctuations [HI, [33[ . Then we match a 2 and b 2 
to reproduce the values (x 2 ) and (y 2 ) from the GLISSANDO 
profiles. Thus, by construction, the spatial rms radii of 
the initial condition and its eccentricity is the same as 
from the Glauber calculation. Nevertheless, the shape 
is not the same, as is evident from Fig. [TJ The Gaus- 
sian profiles are sharper near the origin, which results 
in a faster buildup of the Hubble-like flow in the hy- 
drodynamical stage. Admittedly, the initial density and 
flow profiles should eventually be obtained from the early 
dynamics, such as the Color Glass Condensate theory 
[34 [35L In practice, however, modeling of the par- 
tonic stage carries uncertainty in its parameters. In addi- 
tion, other effects in the early dynamics are present, see 
e.g. [13, [HI], making precise profile calculations difficult. 
Thus the use of a simple parametrization of the initial 
profile is a profitable and practical approach, while the 
need remains for its detailed dynamical justification. 

The Glauber calculations, needed to obtain the a 
and b parameters, correspond to the mixed model [3^ |. 
where the number pf produced particles is proportional 
to (1 — a)N w /2 + aNbin, with N w and iVbin denoting the 
number of wounded nucleons [i(| and binary collisions, 
respectively. The parameter a = 0.145 for RHIC [41, 42] 
and is set by us to 0.2 for LHC. The inelastic nucleon 
cross section is 42 mb for RHIC and 63 mb for LHC [H| . 
The simulations incorporate the fluctuations of orienta- 
tion of the fireball (the variable-axes geometry), which 
result in increased eccentricity compared to the fixed- 
axes geometry [iij . Finally, the expulsion distance of 
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TABLE I: Shape parameters a and b of Eq. (fl) for various centrality classes obtained by matching (x 2 ) and (y 2 ) to GLISSANDO 
simulations, the variable-axes eccentricity e*, and the chosen central temperature T,:. 



c[%] 


0- 5 5-10 10-20 20-30 30-40 40-50 50-60 60-70 70-80 


a [fm] 
b [fm] 

e* 

Ti [MeV] 


RHIC 

2.70 2.54 2.38 2.00 1.77 1.58 1.40 1.22 1.04 
2.93 2.85 2.74 2.59 2.45 2.31 2.16 2.02 1.85 

a r\ o n i o a i o a o c n 01 a a /i 1 c\ A a a c n 

0.08 0.12 0.18 0.25 0.31 0.36 0.41 0.4b 0.52 
500 491 476 460 429 390 344 303 261 


a [fm] 
b [fm] 

e* 

T s [MeV] 


LHC 

2.65 2.47 2.22 1.95 1.73 1.56 1.40 1.24 1.06 
2.89 2.80 2.69 2.55 2.41 2.28 2.15 2.02 1.85 
0.09 0.12 0.19 0.26 0.32 0.36 0.40 0.45 0.51 
768 768 768 768 768 768 768 768 768 



0.4 fm is used in the generation of the nuclear distri- 
butions, and the source-dispersion parameter of 0.7 fm is 
used. This parameter describes the random displacement 
of the source from the center of the wounded nucleon or 
the binary-collision position [3l| . The values of the a and 
b parameters for various centralities and the correspond- 
ing eccentricity parameters 



are collected in Table fl] 

The energy-density profile ([1]) determines the initial 
temperature profile via the equation of state 29]. The 
initial central temperature, Xi, is a model parameter 
dependent on centrality. For RHIC calculations it is 
adjusted to reproduce the total particle multiplicity. 
For the LHC we use for simplicity a common value of 
768 MeV at t = 0.25 fm which corresponds to 500 MeV 
at t = 1 fm. When the LHC multiplicity data are avail- 
able in the future, this parameter will be tuned more 
realistically. 



B. Hydrodynamics 

The hydrodynamic equations used in this work were 
described in detail in Refs. [H, l29| - We use inviscid 
(ideal-fluid) , baryon-free, boost-invariant hydrodynam- 
ics. The equations are written in terms of the veloc- 
ity of sound, c s , whose temperature dependence encodes 
the full information on the equation of state of the sys- 
tem. Importantly, we incorporate the known features of 
c s (T), which at high temperatures are given by the lat- 
tice QCD calculations, at low T follow from the hadron 
gas including all resonances, while in between an inter- 
polation is used. Importantly, no phase transition, but a 
smooth cross-over is built in, in accordance to the present 
knowledge of the thermodynamics of QCD at zero baryon 
chemical potential. The plot of the resulting c s (T) can 
be found in Ref . [II] . 



The initial proper time for the start of hydrodynamics 
is fixed to have the value 

t = 0.25 fm (3) 

for all centralities, both for RHIC and LHC. This early 
start of hydrodynamics allows for a fast generation of 
transverse flow. 1 



C. Freeze-out 

The hydrodynamic evolution proceeds until the freeze- 
out occurs, where the assumed condition for the univer- 
sal freeze-out temperature is T/ = 145 MeV. This value 
is somewhat lower than in several fits of the chemical 
freeze-out [ill [46|, [47|, however, it agrees with the re- 
cently made global fits to particle transverse momentum 
spectra of Ref. [H, H§] , where the value around 145 MeV 
was obtained for the kinetic freeze-out temperature. 2 

The freeze-out hypersurfaces are compared in Fig. [5] 
We use the centrality 20-30% as an example, as for other 
cases the results are qualitatively similar. The left and 
middle panels show the role of the initial condition used 
for the RHIC analysis. For the Glauber initial condi- 
tion used in Ref. (27| (left) the freeze-out hypersurface 
is smaller in transverse size but longer in the life-time 
than for the case of the Gaussian initial condition (mid- 
dle) used in the present work. Parameters in both cal- 
culations were optimized to reproduce the yields, the pt 
spectra, and i»2- The fact that for the Gaussian initial 



The ignition of hydrodynamics may be delayed to later times, 
about 1 fm, if it is preceded with particle free-streaming starting 
at 0.25 fm followed by the Landau matching [24ll , This mech- 
anism generates some initial flow (transverse and elliptic) for 
hydrodynamics. The resulting freeze-out hypersurfaces are very 
close to those used in this work. 
2 The use of this lower freeze-out temperature needs the intro- 
duction of the strangeness inequilibrium factors 7 S in order to 
reproduce the abundances of strange particles |48|I . 
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FIG. 2: (Color online) Sections of the freeze-out hyp ersurf aces for the Glauber initial condition used in Ref. [27| (left) and 
for the Gaussian initial condition (middle) for centrality 20-30% for RHIC, and for the Gaussian initial condition for centrality 
20-30% for LHC (right). The dots with numbers indicate the values of the transverse velocity at freeze-out. 




profile the source has a significantly larger transverse size 
and shorter emission time is a crucial feature for the suc- 
cessful description of the HBT data. 

The freeze-out hypersurfaces have the volume emission 
parts (with time-like normal vectors, i.e. running approx- 
imately flat in Fig. [2]), which arc similar to the blast- wave 
parametrization [50(. The lifetime is short: about 9 fm 
for central (c = 0-5%) and 7 fm for mid-peripheral (c = 
20-30%) collisions at RHIC. However, they also contain 
the surface emission parts (with space-like normal vec- 
tors, i.e. running more-or-less vertical in Fig. [2|), absent 
in the usual blast-wave parametrizations. Nevertheless, 
there is no problem with the non-causal surface emis- 
sion [El], H2, Ell , which is negligible, since less than 0.5% 
of particles are emitted back into the hydrodynamic re- 
gion. This very small fraction follows from the large 
transverse flow velocity at large radii, as indicated by the 
labels in Fig. [2] This large flow carries the particles out- 
ward. In addition, for the RHIC case the hypersurfaces 
are not bent back at low times, as occurs in some hy- 
drodynamic calculations. Quantitatively, we have found 
that about half of the produced particles comes from the 
volume emission part and about half from the surface 
emission part. The relevance of the surface emission is 
also stressed in Refs. [55|. 

The right panel in Fig. [2] shows the freeze-out hyper- 
surface for the LHC extrapolation for the Gaussian initial 
profile and the same centrality. We note a larger trans- 
verse size and emission times compared to the RHIC case 
of the middle panel. The surface-emission parts of the 
curves are bent back, but at the same time the expan- 
sion velocity is larger, such that the backward emission 
problem is again irrelevant. 

THERMINAT0R is used to carry out the statistical 
hadronization at the freeze-out hypersurface according 
to the Cooper- Frye formulation [5g|. According to the 
assumed single-freeze-out approximation, identifying the 



kinetic and chemical freeze-out temperatures, rescatter- 
ing process after freeze-out are neglected. We have 
checked that the collision rate after freeze-out is mod- 
erate for the hypersurfaces applied in this work. We es- 
timate it by considering a pion straight-line trajectory 
and counting the number of encounters with other parti- 
cles closer than the distance corresponding to the pion- 
hadron cross section. The average number of these tra- 
jectory crossings is about 1.5-1.7 per pion. This shows 
that the single-freeze-out approximation |46| works rea- 
sonably well for the present case. At a more detailed 
level, one could use hadronic afterburners to model the 
elastic collisions [13, [H, l59|. or attempt the hydro-kinetic 
approach implemented in [60j] . 

D. Model parameters 

For the reader's convenience, we list here again all the 
model parameters. The model has the total of 5 parame- 
ters. Parameters dependent on the centrality classes are 
the Gaussian widths, a, b, and the initial central tem- 
perature, T t . The widths are fixed with GLISSAND0 (HI, 
while Ti is adjusted in order to reproduce the particle 
multiplicities. The parameters independent of centrality 
are the starting proper time of hydrodynamics, To, and 
the universal freeze-out temperature, Tf. 

E. Two-particle femtoscopy 

The method of dealing with the HBT quantities has 
been thoroughly described in Ref. [3(J, where the reader 
is referred to for details. Here we only note that our 
technique follows exactly the experimental procedure of 
extracting the femtoscopic quantities, including the pas- 
sage to the LCMS (local co-moving system) frame, the 
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application of the two-particle method, and carrying out 
the Bowler-Sinyukov procedure [6l|, [62j to incorporate 
the Coulomb corrections. 

We use the standard Pratt-Bertsch parametrization 
[64l [HI . For the case without azimuthal symmetry the 
correlation function is fitted to the formula including the 
cross terms between the HBT radii 661: 



C(<Zout, 9side, 9long) = 1 + Acxp(-i?^ 

_r,2 2 
-^longyiong 



□2 2 
^sidcysido 



-*^out — sidc^out 9sidc -^out — long Qout Q'loiig 
-^sidc — longQsidcQ'long) (4) 

The basic quantities of the azimuthally sensitive HBT 
analysis are the second-order Fourier coefficients defined 



Rl, tn + 2Rl tto ms(2(f>), 



R al 



side 



(4>) 

-Rout-side (</>) 
-Rout -long 
-Rsidc-long(</>) 



R=i 



out,0 ' ^-"*out,2 ' 
2 



sido.O 
^long,0 

+ 2#L g , 2 cos(2^ 
2i?out-sidc,2sin(2^), 

2i?out-long,2COs(20), 
2i?side-long,2 sin(2^). 

III. RESULTS 
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A. RHIC 

In Fig. [3] we present the summary of our results com- 
pared to the experimental data from the STAR AuAu 
collisions at ^/snn — 200 AGeV [631 ] . The figure shows 
the centrality and pair transverse momentum (fcr) de- 
pendence. For each centrality, associated here with the 
number of participants N pliIt on the horizontal axis, we 
plot the experimental points (filled dots) and the model 
results (empty symbols) . The points from top to bottom 
correspond to kr contained in the bins of 0.15-0.25 GeV, 
0.25-0.35 GeV, and 0.35-0.6 GeV. The top panels show 
the radii squared averaged over the <f> angle, from left 
to right, i^ uti0 , #side,o: and -ftkmg.o- The bottom panels 
show the magnitude of the allowed oscillations divided 
by i?g ido q, which is the adopted convention used in pre- 
senting the experimental data. 

The values of the model points in the plots were ob- 
tained by first solving the hydrodynamic equations and 
running THERMINAT0R as described in the previous sec- 
tions. Then, according to our procedure, 144 separate 
correlation functions have been constructed (6 central- 
ity bins x 4 kr bins x 6 bins for the phi angle). The 
(f> dependence of the correlation function was obtained 
in the following way: an angle between the fc-r direc- 
tion and the reaction plane was determined and the pairs 
were grouped into six bins: (-^, ^), f§), (ff , ff ), 

(if' if)' (if' if)' and (if'ix)- Each of the six size 
parameters was then plotted vs. the reaction plane ori- 
entation for a given centrality and kx bin. The depen- 
dence was fit with the formula ([5]) including only the 



allowed oscillations, and the fit values were corrected for 
the finite bin size [66j. By symmetry arguments, the os- 
cillations in the long, side-long and out-long components 
must vanish, as they do within statistical error in our 
calculations. This is the reason why they are not shown 
in the plot. 

Fig. [3] shows a very good agreement for i? O ut,0 an d 
-Rside.O) where the model points are close the the experi- 
ment at all centralities and kx bins. The agreement for 
i?i 2 ong o is somewhat worse in the lowest kx bin, with the 
model overshooting the data by up to 10% for the val- 
ues of the radius itself. This may be due to the assumed 
exact boost invariance in the model calculation. The os- 
cillations -Rg idc 2 an d ^?out 2 are m a very good agreement 
with the data. The oscillations in the cross term i? ou t-sidc 
are about 50% above the data. This mismatch requires 
explanation. 



B. Gaussian vs. Glauber initial condition 

In order to point out the relevance of the initial con- 
dition mentioned in Ref. (24j , we compare in Fig. [5] the 
angle-averaged side, out, and long radii, and the ratio 
^out/-Rsidc, obtained with the Glauber and Gaussian ini- 
tial conditions. We note a much better agreement with 
the presently applied Gaussian initial condition. In par- 
ticular, the ratio i?out/-Rsidc coincides now with the ex- 
periment, while with the Glauber profile it was about 
20% above the data and the fc^-dependence was not re- 
produced. 



C. Predictions for LHC 

In Fig. [5] we show a plot analogous to Fig. [3] but with 
calculations done for the LHC energies, assuming the ini- 
tial temperature of 768 MeV. As expected, the larger 
radii reflect the growth of the overall size of the system at 
freeze-out. All three radii grow, but the out radius seems 
to grow less, which is discussed in detail later. Also the 
azimuthal oscillations relative to -Rg idc are smaller com- 
pared to RHIC. Again, this is expected: the initial asym- 
metry of the system at a given centrality is quite similar 
at RHIC and at the LHC, since it is mainly driven by the 
overlap geometry of the two colliding nuclei. The sys- 
tem then evolves from the out-of-plane extended source 
towards the more spherical shape, and if the evolution 
time is long enough, it may eventually overshoot and be- 
come in-plane extended. Measurements at RHIC show 
that the source freezes out while still in the out-of-plane 
shape. Our calculations for LHC show that the evolution 
time to freeze-out is longer, but not long enough to pro- 
duce the overshoot. The smaller oscillations (relative to 
i?g ido ) are a consequence of the fact that the system has 
been evolving for a longer time and effectively becomes 
more spherical than at RHIC. 
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FIG. 3: (Color online) Results for the RHIC HBT radii and their azimuthal oscillations. For each value of N pax t on the 
horizontal axis we plot the experimental points (filled symbols) and the model results (empty symbols). The points from 
top to bottom at each plot correspond to kr contained in the bins 0.15-0.25 GeV (circles), 0.25-0.35 GeV (squares), and 
0.35-0.6 GeV(triangles). The top panels show i?o U t,o, ^sidc.oi an d A 2 ong.o> the bottom panels the magnitude of the allowed 
oscillations divided conventionally by i?sid c ,o- Data from Ref. [63| . 
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FIG. 4: (Color online) Results for the angle-averaged side, out, and long radii, and the ratio Rout/ R B ide, for the Glauber initial 
condition (left) and for the Gaussian initial condition (right). Data from Ref. [r33j |. 



In Fig. Owe show an example comparison of the de- 
pendence of the radii at RHIC and LHC energies at one of 
the centralities (in this case top 5%). Two main features 
are apparent: the overall increase of the system size as 
well as larger gradient (steeper slope) of transverse radii 
considered as functions of kr- The latter is a consequence 
of the larger averaged transverse flow developed at the 
LHC. This is illustrated in Fig. [2]- compare the numbers 
indicating transverse velocity at freeze-out on panels b) 
and c). The cause for the former was already discussed, 



now we concentrate on one particular feature. While the 
overall size of the system is indeed larger at the LHC, the 
radii do not seem to grow in the same way. The out radius 
grows significantly less that the other two, which is best 
illustrated by the -R ut/-Rsidc ratio, which decreases from 
1.1 at RHIC to 0.95 at the LHC. It comes from a qualita- 
tive change in the results of the hydrodynamic calculation 
seen in Fig. [2j At RHIC the hypersurface at freeze-out 
temperatures was "outside-in" , or in other words parti- 
cles at the larger transverse distances froze out earlier. 
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FIG. 5: (Color online) Predictions for the LHC. Filled symbols are RHIC data from STAR, shown for comparison (circles - kr 
in 0.15-0.25 GeV, squares - kr in 0.25-0.25 GeV, triangles - kr in 0.35-0.6 GeV). Open symbols are the results of our calculation 
for the LHC energy with the initial temperature set to 768 MeV. The last kr bin has been divided into two: open triangles 
- kr in 0.35-0.45 GeV and open diamonds - kr in 0.45-0.6 GeV. Top panels show sizes (squared) averaged over the (f> angle. 
Bottom panels show magnitude of the allowed oscillations, divided by R^dc- Data from Ref. [Hj]. 
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FIG. 6: (Color online) Comparison of calculations for the top 5% central AuAu collisions at RHIC (squares) and PbPb collisions 
at LHC (circles) energies. Lines are drawn to guide the eye. Top-left panel shows out radius, top-right - side radius, bottom 
left - long radius, bottom right - out over side radii ratio. 



We refer the reader to our previous work [3(| where this modified Blast-wave parametrization. The RHIC data 
effect was studied in the simplified form by analyzing the were found to be consistent with the so-called "negative 
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FIG. 7: (Color online) Calculations for the LHC with the Gaussian initial conditions for 6 centralities (filled circles - 0-5%, 
filled squares - 5-10%, filled down triangles 10-20%, open circles - 20-30%, open squares - 30-40%, open up triangles - 40-50%). 
Upper plots show the ^-averaged radii (from left to right: R ou t, -Rside, R\ong, lower plots show the magnitude of the allowed 
oscillations (from left to right -Rout, 2, ^sidc,2> -Rout-side, 2), all as functions of the pair transverse-momentum kr- 



a" scenario, that is the "outside-in" freeze-out, which is 
also consistent with our more detailed hydrodynamic cal- 
culation, and other hydrodynamic calculations. On the 
other hand the so-called "positive a" , or "inside-out" sce- 
nario was inconsistent with the data mainly due to too 
small i?out- The detailed hydrodynamics calculation for 
LHC energies shown in this work exhibit a qualitative 
change in the freeze-out shape to the "inside-out" type. 
Therefore we do expect to see a smaller, or more pre- 
cisely, less increased outwards radii, and that is exactly 
what we see in Fig. [51 The effect is large, when compared 
with the expected experimental systematic uncertainties 
and therefore can be easily tested. Its confirmation in the 
data would be a very strong indication of the existence 
of the hydrodynamically behaving medium in relativis- 
tic heavy-ion collisions, and would indicate that we are 
significantly advanced in understanding the dynamics of 
this medium. 



D. LHC results as function of kr 

In Fig. [7| we show the azimuthally sensitive HBT re- 
sults for the LHC energies as a function of the pair trans- 
verse momentum kx, including the low momentum bin 
of 0.05 — 0.15 GeV. The out and side radii fall with kx 



in a similar way, linearly. The expected decrease of the 
size with growing centrality is also apparent, while the 
Rout /Rside ratio is close to a constant for all centralities 
and values of kx- The long radius also falls as expected. 
In general, the trends in the radii, plotted both versus 
centrality and kx, are smooth and well understood. It is 
interesting to note the behavior of the oscillations. Their 
sign is the same as for LHC as for RHIC. Although we 
do see a decrease in the magnitude of oscillations with 
lowering kx, especially in the side direction, they do not 
change sign. Also the magnitude of the oscillations grows 
consistently with increasing centrality for all kx bins. 



The oscillations obtained from the analysis of the HBT 
correlation functions correspond to the system asymme- 
try at the final stages of the collision, when the freeze-out 
occurs. It is interesting to compare them to the initial 
space asymmetry, as obtained from the Glauber calcula- 
tions and presented in Tab. HI This is shown in Fig. [5] 
One can see that after dividing the observed asymmetry 
by its initial value one obtains a curve which appears to 
be universal, within the statistical uncertainties of this 
study, for all considered centralities. 
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(Color online) Scaled azHBT oscillations for the LHC initial conditions for 6 centralities (filled circles - 0-5%, filled 
- 5-10%, down triangles 10-20%, open circles - 20-30%, open squares - 30-40%, up triangles - 40-50%). The symbols 



FIG. 8: 

squares 

show the magnitude of the allowed oscillations (from left to right: R% 
e from Table U and by R^ id 



out,2i -^sidc.25 -ftout — side, 2 



divided by the initial eccentricity 
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FIG. 9: Evolution of the shape of the system with time for the LHC energy and centrality 20-30%, as marked with the produced 
primordial pions. Each panel shows the birth places of primordial charged pions emitted in a time interval of duration 1 fm. 
The pions originate from the region placed at the center of the collision (\z\ < 1 fm). The top-left panel corresponds to the time 
interval 1.0-1.7 fm, top middle-left to 1.7-2.7 fm, and so on up to the time interval 10.7-11.7 fm shown in the lower-right panel. 
The shades of gray indicate the relative number of primordial charged pions emitted from the area element in the transverse 
(x, y) plane within the given time interval. 



E. Emission history the particles are emitted during the whole evolution time, 

the observed shape is essentially a multiplicity-weighted 

The above considerations constitute an overall picture 
that can be deduced from the radii alone. However, since 
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FIG. 10: Same as Fig. [9] for all pions, primordial and from resonance decays. 



average of the shape evolution of the system at subse- 
quent times. An example of such shape evolution with 
indicated birth points of the primary pions is shown in 
Fig. M for the 20-30% centrality for the PbPb collisions 
at LHC. The information discussed below is a more de- 
tailed analysis of the behavior which may be inferred 
from Fig. [2] The panels in Fig.[9]show where the charged 
primordial pions are emitted in a given time interval. The 
scale (shown in the bottom-right panel) indicates how 
many pions are emitted from a given area element in the 
transverse plane. The picture shows the slice centered 
at midrapidity (the z coordinate of the particle emission 
satisfies the condition \z\ < 1 fm). At the very beginning 
of the collision (t < 1.0 fm) there are no particles emitted. 
However, we observe particle emission at the very early 
stage - already in the first time slice (1.0-1.7 fm). Obvi- 
ously the system has only begun its evolution, hence the 
shape is still very much out-of-plane extended, reflecting 
the initial overlap geometry. We also see that up to the 
very last evolution stage the emission is only from the 
surface of the system. As the evolution progresses par- 
ticles are emitted from more and more spherical shells, 
up to the time of 9.7 fm. At this moment system reaches 



the freeze-out temperature in its full volume and emis- 
sion changes its nature from the surface-like shell emis- 
sion to the bulk emission from the whole volume. It is 
especially interesting to look at panel 10 (9.7-10.7 fm). 
One can see that the volume emission starts in the cen- 
ter of the system first, this is a very good illustration of 
the "inside-out" type of emission discussed earlier. What 
is also interesting is that even though the overall shape 
of the system is very much in-plane extended, the in- 
ner emitting part has an out-of-plane shape. In the final 
stage of the evolution (panel 11: 10.7 - 11.7 fm) the vol- 
ume emission continues and finally the outer parts of the 
system emit particles. As can be seen from the picture 
and the discussion, the emission patterns coupled to the 
shape evolution of the system are quite complex. There 
is no way to tell, a priori, what will be the overall system 
shape seen by femtoscopy. One must perform a detailed 
simulation of all the evolution stages, which must also 
include a realistic simulation of the number of particles 
coming from each stage of the evolution. This is precisely 
what we have done in this work. 

In Fig. [10] we show the analog of Fig. [9l but with all 
pions, not only the primordial ones. We note that the res- 
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onance decays "wash out" the production regions, which 
is effectively increasing the size. This is due to the trans- 
verse flow, which has a tendency to carry out the reso- 
nance outward before if decays into pions. Qualitative 
conclusions concerning the surface and volume emission 
are the same as for Fig. GO We also note that pion emis- 
sion continues well past the the time of 11.7 fm, as reso- 
nances with longer lifetimes subsequently decay. The de- 
tailed discussion of the influence of strongly decaying res- 
onances on the HBT correlation functions and extracted 
radii can be found in (30j . 

IV. CONCLUSIONS 

The results presented in this paper show that it is pos- 
sible to achieve a uniform description of the RHIC soft- 
hadronic data including azimuthally sensitive HBT radii. 
The initial space asymmetry of the source, tuned earlier 
to reproduce the value of t>2, turns out to be precisely 
such that the azimuthal dependence of Rg idc 2 and i?o Ut 2 
is also very well described. This verifies good consistency 
of our approach consisting of relativistic hydrodynamics 
and statistical hadronization - only the model prediction 
for the out-side radius is significantly different from the 
data, which requires explanation. For the LHC energy 
we predict the moderate increase of the HBT radii and 
the decrease of their azimuthal oscillations. Such effects 
can be naturally explained by longer evolution times at 
LHC. The space-time patterns of particle emission were 
discussed in detail. They indicate that the shape of the 
system seen by the femtoscopic methods is an average of 
the complex and varying in time shape of the emitting 
source. 

In summary, the results of our calculations for RHIC 
and LHC conditions show several notable features. First 
of all the calculations for RHIC show a remarkable agree- 
ment with the broad spectrum of the soft physics data. 
We have concentrated on the transverse dynamics of the 
source, and our model appears to properly describe not 
only the momentum part of the observed phase space, but 
also the space-time part. The unique feature of the model 
is the proper description of the azimuthal asymmetry in 
the side and out directions, again both in the momen- 
tum and space-time, which has been achieved before only 
in simplified and non-dynamical blast-wave parametriza- 
tions [231 ] which neglect the important contributions of 
the strongly decaying resonances, as well as the surface 



emission. The underlying hypothesis of our work is that 
the system created in heavy-ion collisions at RHIC be- 
haves as a single piece of matter, and can be described 
by the hydrodynamic equations. These equations include 
the state-of-the-art equation of state which assumes that 
the matter above a certain critical temperature is in a de- 
confined phase (Quark Gluon Plasma). In our calculation 
the evolving system at RHIC spends a significant amount 
of time in that phase, and its properties are essential in 
shaping up the final observables. This work provides an- 
other crucial confirmation that such a hypothesis is con- 
sistent with the experimental data. While in itself it does 
not constitute a proof that the Quark Gluon Plasma is 
indeed created at RHIC, we stress that any alternative 
explanation must at least achieve a similar agreement 
with the experimental data to be considered viable. 

The system created in heavy ion collisions at the LHC 
is predicted to be in many ways similar to the one cre- 
ated at RHIC, at least in the sense that it also is ex- 
pected to spend a long time in the deconfined phase. In 
fact, this time is predicted to be significantly larger than 
at RHIC, such that its influence on the final observables 
may be more pronounced. Again, we work under the as- 
sumption that hydrodynamics provides a good descrip- 
tion of the matter in such conditions, which allows us 
to provide predictions for final state observables at these 
energies. We stress that the underlying mechanisms of 
the model do not change at all between the RHIC and 
the LHC energies, only a few of the external parame- 
ters, such as the initial nucleon-nucleon cross-section or 
the initial temperature are changed in a reasonable way. 
Nevertheless, we are able to identify significant changes 
in the observables (with respect to RHIC) that can be 
easily measures in the LHC experiments. In particular 
in this work, which focuses on femtoscopy, we have iden- 
tified two of them which are particularly sensitive: the 
decrease of final observed anisotropy of the source (rela- 
tive to i?g idc ), reflected in the decrease of the oscillations 
in R ou t and R S ide radii, and the change from the "inside- 
out" to the "outside- in" type of freeze-out, reflected in 
the femtoscopic radii themselves and best illustrated as 
the decrease of the R ou t / 1 Rside ratio. These features are 
specific enough such that they provide strict tests of the 
validity of the hydrodynamic hypothesis. If observed, 
they will be a strong argument that systems at RHIC 
and the LHC can indeed, at least in the soft sector, be 
described by essentially the same physics principles. 
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